Evolution of the quasiparticle spectral function in cuprates 
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| | ' We analyzed photoemssion data for several doping levels of the Bi2Sr2CaCu20$+ x compounds, 

ranging from overdoped to underdoped. We show that the high frequency part of the spectra near 
(0, 7r) can be described by Fermi liquid theory in the overdoped regime, but exhibits a non-Fermi 
^ , liquid behavior in the underdoped regime. We further demonstrate that this novel behavior fits 

S_h ■ reasonably well to a \j \[lo behavior suggested for systems with strong spin fluctuations. 

^ ' 74.20.Jb 



One of the most intriguing characteristics of the cuprates is their qualitative change in electronic properties with 
doping. These changes have been attributed mostly to the pseudogap that appears in the normal state at frequencies 
comparable to the superconducting gap. The goal of this communication is to show that in addition to the pseudogap, 
the high frequency part of photoemission spectra exhibit qualitative changes with decreasing carrier concentration. 
We argue below that this evolution with doping indicates that the cuprates are near a quantum critical point and 
O , exhibit a non-Fermi liquid spectral function for high frequencies. 

A series of recent papers have reported that the spectral weight of the quasiparticle peak near the Fermi surface 
becomes smaller for underdoped samples, and that there is a transfer of spectral weight to the high-frequency part 
of the spectrum |J. In this paper, we present a functional analysis of such data. Our main result is that the 
photoemission intensity I(u>) undergoes a particular type of crossover behavior with decreasing doping (since energy 
is hu>, we subsequently refer to the frequency to). The spectral function changes from a behavior which can be 
reasonably well described by the Fermi liquid form J -1 (a;) oc (r 2 + to 2 ) to a behavior which to good accuracy can 
be described by a non- Fermi liquid form / _1 (cj) oc ^/lo. Such behavior is predicted when the dominant effect of 
quasiparticle decay is the interaction with overdamped spin fluctuations Q . 

We begin by reviewing the theoretical aspects of the problem. In the sudden approximation, the photoemission 
intensity at a given momentum k is given by Ik(Lo) = Ak(j-o)n(Lo) where Ak(t-u) = (l/ir)ImGk{<jj) is the quasiparticle 
spectral function, and 12(10) is the Fermi function which selects states occupied by electrons. Fermi liquid theory 
predicts that, regardless of how strong the interaction is, the single-particle Green's function at low frequencies has 
I ■ the form 

c ' z 

(lo - (e k - e F ) + iu)\u)\/T) 

where tp is the Fermi energy. For momenta near the Fermi surace, i.e., at small — ep, the spectral function which 
emerges from Gu(lo) possesses a quasiparticle peak at lo = — tp\ the width of the peak decreases as (e^ — £i?) 2 
rS \ with the approach of the Fermi surface. At the Fermi-surface, the quasiparticle peak transforms into a ^-function for 
lo — 0, while for frequencies lo < T, Ak{uj) takes the form A^. 1 {lo) oc {lo 2 + Y 2 ). The inverse proportionality to to 2 
at low frequencies is a fundamental property of Fermi liquids which does not depend on details of the quasiparticle 
interaction. However, T, which serves as the upper frequency cutoff for this universal behavior, is model dependent. 
When the Fermi system is far from any instability, T is of the order as €f, which is typically the upper cutoff for a 
low energy description (i.e., expansion in powers of lo). This changes if the system is close to a quantum critical point 
where it undergoes a spontaneous symmetry breaking. In this situation, T <C cf and T scales with the distance from 
the critical point. For uj <§; T, the system still possesses Fermi liquid behavior- i.e., A^ 1 ^) oc (lo 2 + T 2 ). However, 
for r <^ lo <^ tp the system crosses over into the region which is in the basin of attraction of the quantum critical 
point. In this region the system behavior is governed by critical fluctuations. These fluctuations impose a novel 
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frequency dependence of the spectral function up to frequencies comparable to the Fermi energy. It is essential that 
this picture survives even if the system trajectory in the phase space of parameters slightly bypasses the quantum 
critical point (e.g., the actual transition into the symmery breaking state can be a weak first order transition or 
involve an intermediate phase) or the system possesses a low-energy crossover from one type of critical behavior to 
that associated with small but relevant perturbations around the quantum critical point (in cuprates this may be 
due to a Fermi surface evolution very near the magnetic instability). These effects reduce the range of Fermi-liquid 
behavior to even smaller frequencies, but they do not modify the behavior at u) > T as long as V remains the largest 
energy scale associated with instability. 

In recent years there have been a number of suggestions for possible quantum critical points in the cuprates Q . In our 
opinion, the most plausable candidate is the instability point towards antiferromagnetism. The system approaches 
this critical point as the carrier concentration decreases, and the system is antiferromagnetic at half-filling. The 
quantitative measure of closeness to a magnetic transition is the spin correlation length, £. When £ is large enough, 
the system behavior is predominantly critical over a wide frequency range. As £ clearly increases with decreasing 
doping, the range of critical behavior also increases. 

Recently, one of us (AC) considered the interaction between fermions and overdamped spin fluctuations and found 
that r ~ oj s f where oj s f oc £~ 2 is a typical spin-fluctuation scale 0]. More specifically, the fermionic Green's function 
along some portion of the Fermi surface near (0, tt) and related points was shown to have the form 

G -i( w ) oc 1 ^ (2) 

z 1 + v 1 - i\uj\/uj a f 

where Z oc At small frequencies, u> < ui s f, this Green's function has a conventional Fermi-liquid form G~ 1 (lu) oc 
(u) + iuj\uj\/(4u) s f)). At higher frequencies, however, the Fermi-liquid form crosses over to a region of novel frequency 
dependence- G(u>) oc Ae~ 47r / 4 |aj|~ 1//2 s<?n(jj, where A — Zlo s j^ 2 . As A is independent of £, it remains finite when 
£ — > oo, as it should for quantum critical behavior. One can probe this behavior at finite £ for \u\ Values of 

u) s f for various dopings have been extracted from available NMR data Q|. The NMR data indicate that ui s f (at a 
given T) increases with doping and remains ~ 20 — 30meV (i.e., small compared to cf) for overdoped cuprates. From 
this perspective, some features of the crossover to ^Juj behavior should be observed even in the data for overdoped 
cuprates. 

The tendency towards quantum critical behavior is much less pronounced near the Fermi surface crossing along the 
zone diagonal ((7r, tt) direction). Here, calculations show that the analog of uj s f remains finite for £ — > co with a value 
comparable to sf- The system is therefore in the moderate coupling regime; A(ui) at intermediate u> should display 
a frequency dependence which is intermediate beween quantum-critical and Fermi liquid forms. 

We now turn to the analysis of the photoemission data. The photoemission data were measured using angle-resolved 
photoemission equipment at the Wisconsin Synchrotron Radiation Center and at Ecole Polytechnique Federate Lau- 
sanne; the equipment details are provided in Ref. Samples were fabricated as reported in the literature H, and 
were transferred via a load-lock system so as to retain the oxygen content. The data were taken at a temperature 
well into the normal state (typically 2QQK for underdoped samples, 1QQK for overdoped samples), and above the 
temperature, J 1 *, below which a pseudogap would be observed. In addition, we analyzed the data from Ref. ffl] and 
found our conclusions to be consistent with that report. Altogether, we analyzed data on seven doping levels, includ- 
ing overdoped (T c « 52K, 65K and 75K), optimally doped (T c « 90K), and underdoped (T c « 60K and 30K), in 
addition to analyzing underdoped data (Tc ~ 85K) from Ref. [Q. The details of all samples and doping levels are 
reported elsewhere M . 

We analyzed the photoemission data using seven different techniques, including fitting using the chi-squared criterion 
(with four different background subtraction methods), the Kolmogrov-Smirnov criterion, the ratio of data at different 
binding energies, and the inverse of the data. We provide a complete, detailed discussion of these analysis methods 
elsewhere All seven methods led to the same conclusion, that there is a qualitative change in the nature of the 
spectral function with doping. Due to space limitations, we emphasize two of these methods in this report. 

We first summarize our results. For overdoped cuprates near (0,7r), we found that the data below ~ 250meV 
(the exact value depends somewhat on the background subtraction method employed) were consistent with the Fermi 
liquid Green's function (Eq. [l]) and were inconsistent with the Green's function of Eq. |[ For higher frequencies, 
however, we found some evidence for a crossover to quantum-critical behavior. For overdoped cuprates along the 
(71", 7r) direction, the data were also consistent with Eq. [j]. For underdoped samples near (0, tt), we found that the 
data were consistent with Eq. and were inconsistent with a Fermi liquid Green's function over the entire frequency 
range studied. In the (it, tt) direction, the data are fit equally well by either the Fermi liquid or the non-Fermi-liquid 
spectral functions. 
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We now discuss in more detail how we actually fitted the data. We began by checking whether the finite energy 
resolution of our equipment significantly modified the functional forms of either Eqs. |l| or ||. To do this, we generated 
spectra by convoluting each of the two theoretical spectral functions with the experimental resolution function, and 
subsequently fitted the results using the unbroadened functional forms of Eqs. [I] and |[ We found that the finite 
energy resolution of the electron analyzer has only a small impact on the analyses and could be almost fully absorbed 
into small renormalization of the input parameters, ui s f and T. 

In our first method, we took special care in dealing with the background. The key assumption we made is that pho- 
toemission spectra for wavevectors far from the Fermi surface are indistinguishable from the background. Analyzing 
the data in these fc-ranges, we observe that they are nearly frequency independent over a relatively wide frequency 
range. The insets of Fig. [I] illustrate this situation (the flat regions are indicated by arrows). We further assumed 
that the flat background is independent of k, i.e., it remains the same near the Fermi surface. We then eliminate this 
background by replotting the data in the form R(u>) = — I (u m ax)) / (I (w m in) — where uj m i n and u> max are 

the lower and upper boundaries of the region where the background is independent of frequency. Plotting the data 
in this form additionally allows us to eliminate the overall factor in the intensity. Finally, we fit R(lo) by Fermi-liquid 
and spin-fluctuation forms using oj s f and T as adjustable parameters and check which form works better. 

Applying this method to experimental data, we obtain the results presented in Fig. Jl|. The actual spectra from 
which the fits were obtained are the upper curves in the insets of each part of Fig. We include three sets of 
data- oxygen overdoped material (Tc ~ 65K), oxygen underdoped material (Tc ~ 30K and 60K), and Dy-doped 
underdoped material (Tc ~ 85K) from Ref. 0. The spectra analyzed are those for the /c-points closest to the Fermi 
surface. For overdoped samples near the (0, 7r) point, the data consistently yield a better fit for Eq. [j]. Along the 
(tt, tt) direction, the data are also fit better by Eq. [j], although the error bars are more significant. For underdoped 
materials near the (0, tt) point, the spectral function from Eq. |^ yields a much better fit than that of a Fermi liquid. 
Furthermore, we consistently obtain a spin fluctuation energy of 30- 40 meV, which is a bit larger than, but still 
comparable to, the to s f ~ 10 — 20meV extracted from NMR. Along the (tt, tt) direction, both spectral functions fit 
the data equally well, thus we cannot judge which form works better. 

In the second method of analysis, we assumed that the background is independent on k over the whole frequency 
range studied in photoemission. After making this assumption, we subtracted the data for fc-values far away from 
the Fermi surface from the data near the (0,7r) point, inverted the subtracted data, and fitted the results to both 
Fermi liquid (uj 2 ) and \fuj forms. This method allowed us to fit the data over a wider frequency range than the first 
one. The results are shown in Fig. |2|. The results provide additional insights to the results in Fig. [I]. For overdoped 
samples, we find that the data are fit best by a quadratic (Fermi liquid) form for binding energies up to ~ 250meV, 
and by a y/uj form for higher binding energies. The exact binding energy of the crossover depends somewhat on the 
background subtraction method employed. For underdoped samples, the data in Fig. ^| are better fit by a y/U form 
for the entire binding energy range above ~ lOOmeV. We also performed the same analysis for the data taken near 
the Fermi surface crossing along the (tt, tt) direction (not shown in Fig. ^|). For both overdoped and underdoped 
samples, the Fermi liquid and non-Fermi liquid forms fit the data equally well. As a consequence, it is difficult to 
assess which form works better in the along the (tt, tt) direction for any doping level. 

In summary, theory predicts that when doping decreases, the Fermi-liquid form of the spectral function near (0, tt) 
exists in progressively smaller range of low frequencies, while at higher frequencies the spectral function crosses over 
to a novel frequency dependence, A _1 (w) cx y/uj, associated with the closeness to an antiferromagnetic quantum 
critical point. Experimentally, we have measured and analyzed the high frequency part of the photoemission data for 
BSCCO-22I2 materials and its variation with doping. We found that the data near (0,7r) for underdoped cuprates 
agrees well with the spin-fluctuation form of the spectral function, Eq. ||. For the frequency range we used, it yields a 
non- Fermi liquid form j4 _1 (oj) oc v/M". I n contrast, the data for overdoped cuprates is well- fitted by the Fermi liquid 
form A^ 1 ^) cx (lu 2 + T 2 ) over a substantial frequency range (up to ~ 250meV). At larger frequencies, we found some 
evidence for crossover to I/^/ja7[ behavior. Overdoped cuprates in the (tt, tt) direction exhibit Fermi liquid behavior. 
Neither analysis method yields conclusive results for underdoped cuprates in the (tt, tt) direction. 
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FIG. 1. The fits of the normalized photoemission intensity R{lj) by the Fermi-liquid frequency dependence, Eq. [I] (dotted 
line), and by the spin-fluctuation form, Eq. |^ (solid line). The horizontal axis is the electron binding energy from (the Fermi 
energy) and the hash marks 100 meV apart. The vertical axis is R{u) on a logarithmic scale. Figs, (a)-(c) are fits to the data 
near the (0, ir) point including: (a) oxygen underdoped, Tc ~ 60K material, (b) Dy-underdoped, Tc ~ 85ii" material from 
Ref. j?]], and (c) oxygen overdoped, Tc ~ 657^ material. Figs (d)-(f) are fits to the data near the Fermi surface crossing in the 
(n,n) direction, including: (d) oxygen underdoped, Tc ~ 307^ material, (e) Dy- underdoped, Tc ~ 85K from Ref. Q (same as 
in Fig. (b)), and (f) oxygen overdoped, Tc ~ 65K material. Insets to the figures show angle-resolved photoemission spectra 
taken near the Fermi surface (upper spectra) and far away from the Fermi surface (lower spectra). The horizontal axis is the 
electron binding energy, from (Fermi energy) and the hash marks are 100 meV apart. The arrows indicate the frequency 
range for which the spectra far away from the Fermi surface are frequency independent. This is the range over which we fitted 
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FIG. 2. Results of the second analysis method of fitting the data. We substracted the data far away from the Fermi surface 
(which we identify with the background) from the data taken near the Fermi surface crossing near the (0, it) point for all 
fc— values and inverted the difference . We then fitted the inverse intensity by the Fermi-liquid, w 2 form, (dotted line) and the 
^/To form (solid line), which is the high-frequency limit of the spin-fluctuation form (Eq. |2J). Figs (a) and (b) are the fits for 
oxygen overdoped, Tc ~ Q5K material and oxygen underdoped, Tc ~ Q0K material, respectively. Observe that the uj 2 form 
does not fit the data for underdoped sample at all. For overdoped sample, the uj 2 form fits data for binding energy up to 
~ 250meV, above which there is a crossover to the y/u dependence. 
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